Nonharmonic oscillations of nanosized cantilevers due to quantum-size effects 
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Using a one dimensional jellium model and standard beam theory we calculate the spring constant 
of a vibrating nanowire cantilever. By using the asymptotic energy eigenvalues of the standing 
electron waves over the nanometer sized cross section area, the change in the grand canonical 
potential is calculated and hence the force and the spring constant. As the wire is bent more 
electron states fits in its cross section. This has an impact on the spring "constant" which oscillates 
slightly with the bending of the wire. In this way we obtain an amplitude dependent resonance 
frequency of the oscillations that should be detectable. 

PACS numbers: 62.23.Hj, 62.25.-g, 62.30.+d 



I. INTRODUCTION 

Nanoelectromechanical (NEMS) systems are usually 
still in the macroscopic regime in the sense that quan- 
tum effects do not play a large role. This is due to the 
smallness of the quantum energy of a mechanical struc- 
ture even if it is of nanometer dimensions. For example, 
Huj for a cantilever with a large resonance frequency of 1 
GHz will only have an energy of 4 /ieV, which is much 
smaller than the termal energy at normal conditions. The 
corresponding amplitude A of such quantum oscillations, 
given by \kA 2 — Huj, where k is the transverse spring 
constant, will be very low and hard to detect^. How- 
ever, in certain circumstances quantum size effects have 
a significant effect 2 . For example, Stafford et al£ and 
others^—, calculate the tensile force in a nanowire during 
its elongation. They found jumps in the force correspond- 
ing to different numbers of electron states that fits in the 
wires cross section when the wire was stretched. This ef- 
fect has been measured^. The canonical component in 
NEMS is the cantilever, which is used in a large number 
of systems including atomic force microscopes and can- 
tilever based sensors. It is not obvious that bending a 
cantilever will give the same effect as in the Stafford sys- 
tem where the diameter of the nanowire could be reduced 
to a fraction. 

Here, we show that quantum size effects should be in- 
cluded for thin cantilevers, and even if the effect is small, 
at resonance the effect will be detectable due the the high 
accuracy at which frequency of oscillation can be mea- 
sured: down to tenth of mHzi^ We use the same kind of 
free electron model as in earlier studies to calculate the 
spring constant of the cantilever. 



II. 



MODEL 



When a nanowire is bent the length of a fibre on the 
upper side of the wire increases, corresponding to reduc- 
tion of the width of the cross section. Conversely the 
length is reduced on the lower part corresponding to an 
increase of the width of the cross section. In the middle 



is an unaffected neutral line. This effects will change the 
cross sectional area of the nanowire which in turn change 
the tensile force in the wire as it becomes more bent. 
Because the wire is assumed to be thin, only a few wave 
modes under the Fermi level fits in the cross section. We 
use a one dimensional density of states in a similar way 
as in earlier works 3- -. We consider a straight nanowire, 
which we bend by applying a perpendicular force at the 
free end. The other end of the wire is attached to bulk 
material letting electrons flow in and out of the wire. By 
using standard beam theory we find that bending of the 
wire yields an increase in the cross section area. The 
wire has the length L and a quadratic cross section with 
undeformed side do, see fig [TJ When we applying a per- 
pendicular force F at the end of the cantilever at x — L, 
we obtain 



M 



F(L-x), 
M 

T z ' 

o z =0, 



(1) 
(2) 
(3) 



where M is the moment of the force F at position x, I 
is the areal moment of inertia of the cross section and 
, <r z are stresses in the x, y and z direction respec- 



tively. For strains in the x, y and z directions we have 
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where v is the Poisson's ratio and E is the Young's mod- 
ulus. This lead to the deformation of the cantilever cross 
section as shown in fig [TJ For the deformations u v we 
have 



du v 
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vF 
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(L - x)yz + H (x, z) , 



(7) 
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FIG. 1: A cross section S of the cantilever before and after a 
perpendicular force is applied on the wires end. Inset: Side 
view and top view of the cantilever nanowire. The expansions 
and contractions of the cross sections are larger nearer the 
fixed end if the cross section can contract and expand freely. 



using partial integration, where H (x, z) = by symme- 
try. For the deformations u z we have 



du z _ 

dz 
vF 



2EI 



(9) 

(L - x) z 2 + G (x, y) . (10) 



The shear strain j yz is assumed to be zero, so 



du y du z 

7 - = aT+a^ ' 



(ii) 



yielding an expression for G (x, y) to use in equation (10). 
We finally obtain 
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(L-x)yz, (12) 
(L-x){z 2 -y 2 )+f(x), (13) 



where / (x) is the engineering beam theory solution de- 
scribing the bending of the neutral line y = z = which 
only implies a translation of the whole cross section and 
therefore can be left out when calculating its deforma- 
tion. 

When the wire is deformed by the force, the side of the 
cross section is bent inwards with an angle 9/2 on the 
upper part of the wire cross section and is bent outwards 
with the same angle on the lower side. On the top side 
there is a compression of the wire and on the down side 
there is an elongation. The area after the deformation 
is the difference between two circle sectors with radii R 



and r respectively and 9 is the top angle: 
S 
R 



da 

T + A ' 
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(14) 
(15) 
(16) 
(17) 



resulting in S ~ df t (1 + A/£). From (13) we see that 



2EI K ' \ 2 



The radius of curvature I is given by 



^ = — {L-x). 



I 1 dy 2 1 EI 
We then obtain using (18) and (19): 



S = di 



1 / isF(L-x)d 
EI 



(18) 



(19) 



(20) 



A beam of length L with a perpendicular force F applied 
at the end bends a distance Z: 



Z = 



FL 3 kZL 3 



3EI 3EI ' 



(21) 



where F = kZ , and k in turn is the transverse spring 
constant of the wire. We then obtain 
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1 + 



9i^ (L - xy Z 
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(22) 



which is independent of the Young's modulus E. The 
eigenenergies for the standing electron waves that fits 
the cross section is in the limit of large eigenvalues^ 



h 2 4tt 
En = 2^Y n > 



(23) 



where n is the quantum number and m the electron mass. 
The grand canonical potential of the electron gas in the 
nanowire is for low temperature given by&£ 



fl = -^£^ ] ^(E F -E n (x)f 2 dx. (24) 
Integration of (24) using (22) for small bending Z yielded 
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FIG. 2: Example of spring constant k variation as a function 
of bending Z' — \^\ vd ^ obtained from an equation having 
the same type of dependency on n and Z as (32) has, i.e. 
equation (38). The three curves show the forms of the curve 
obtained from (38) for -g^ below, at and above an integer. 
We see that (38) is independent of the sign of Z as is to 
be expected by symmetry. The k-axis always intersects the 
curve at a local maxima or minima. However, this maxima 
or minima may not be very wide depending on how close we 
are that a new electron state will fit when bending the wire 
slightly. We see in the figure that the width is small when 
-g^ is close to an integer (middle curve) and broader farther 
from an integer. 



Using (29) in (24) and taking the derivative yielded 



k = 



where 
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Because (30) should yield the same as (28) for small 
bending, (30) should be corrected with factor 1/3 to ac- 
count for the variation of S with x: 
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(32) 



This expression is valid for arbitrary bending Z . Plotting 
k from equation (32) as a function of Z is shown in fig [21 
Replacing a sum by an integral we found 
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(33) 



Using (28) for small bending and (33), neglecting the 
second term in the RHS of (33) because N > 1, we obtain 
the spring constant for the unbent wire 



which is Blom et air plus a term proportional to the 
down bending Z squared. We have 



E n 

N 



h 2 4tt 

2^" = ^' 

Ep 

Eo' 



(26) 
(27) 



fc 



\f2v 



2 ld 2 m 3 E 5 F 



5tt 2 



(34) 



due to an increase in the electron gas density of states 
that fits in the nanowire when bent. The k-axis always 
intersect the curve at a local maxima or minima of the 
k- values, as we see in fig [2] We may therefore have k = 
fc + const x Z 2 for not too large bending, see (39). 



where N, if rounded off down to an integer, is the number 
of energy levels below Fermi level. Using for example 
do = 4 nanometer and Ep — 5.5 electron volt (gold, 
silver) we obtain the number of energy levels below Fermi 
level N = 183. d Q = 1 nm yields N = 11. The force F 
which is assumed to be due to the electron gas is given 
by F — — ^ = kZ, so we obtain, assuming ^ = (no 
elongation of the wire), the spring constant 



A' 



2m 



2L 3 V TT 2 H 2 



\J Ep — E nn E n( ). (28) 



Disregarding the variation of S with x in (22) and as- 
suming the same deformation in every cross section as at 
the fixed end to simplify the integration we obtain 



S = di 



Qv 2 d 2 Z 2 



8L 4 



(29) 



III. 



DISCUSSION 



We may then rewrite the harmonic equation with a 
bending dependent spring constant: 



d 2 Z o ~ . ^ dZ . . , 

— T + ujIZ + f3Z 3 + 2 7 — = F sin (tot) . 
dt z dt 



(35) 



This is called the Duffing equation. The amplitude A and 
phase c/> of the stationary solution to the linear (j3 = 0) 
differential equation is given by 



.4 
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(36) 



(37) 



Because the weak non-linearity the resonance frequency 
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FIG. 3: The amplitude response of a driven weakly nonlinear 
harmonic oscillator. /3 = corresponds to the harmonic case 
where the resonance frequency is loq. When the amplitude 
is increased (or f3 is increased), the resonance frequency is 
shifted towards higher values for positive /3 and towards lower 

values for negative /3, from cjq to u^ 2 , + jfiA 2 . 



FIG. 4: P/{Cul) as a function of N = ^ obtained by divid- 
ing the second term in (39) with the first term and identifying 

9v 2 d 2 

the spring constant in equation (35). C = g L4 ° . We note 
the singularities at integer numbers, corresponding to when 
we have k — ko + const x\Z\ instead of k = ko + const x Z 2 
which is the normal behaviour near Z — 0. 



is shifted from ujq to \Juq + j/3A 2 , see reference^. We 

then obtain a shift in the amplitude maximum towards a 
higher frequency for positive /3, as shown in figO Because 
frequency can be measured at high precision, even small 
changes can be experimentally detected. The oscillation 
in the weakly non-linear case takes place around the same 
point of equilibrium as in the low amplitude, i.e. linear 
harmonic case. 

We see from fig [2] that j3 is positive for the upper curve 
and negative for the lower curve around Z — 0. The mid- 
dle curve is problematic: very close to Ef/Eq being an 
integer we have k = k a + constx \ Z\. This discontinuity in 
the first derivative of the curve at Z = vanishes however 
quickly when we move away from Ef/Eq = integer. 

The results in this paper should be valid for wires with 
small enough diameter so that the one dimensional distri- 
bution function is a good approximation. From the result 
presented, it should be possible to experimentally deter- 
mine how close we are to that a new state would fit in 
the nanowire. By first measuring the low amplitude fre- 
quency of the oscillating nanowire and then increase the 
amplitude it should be possible to determine if the reso- 
nance frequency is increased or decreased corresponding 
to positive or negative j3. To illustrate the changes in 
equation (32) with different values of N = Ef/Eq we 
made plots. The function to be plotted is 



jv(i+cz 2 ) 

k'{z)= J2 J - x 



1 + CZ 2 (i + cz 2 ) 



2 ■ 



(38) 



where C = 9 g L 4° , which has the same type of dependency 
on n and Z as equation (32) has. We then obtain curves 



as in fig [2l The lowest curve is obtained for N = 5.9, 
the middle curve for N = 6.0 and the upper curve for 
N = 6.1. For small bending Z, (38) becomes 



N 



k'(Z) = ^2 n ^ N ~ n + 

n=l 



(39) 



Dividing (38) with the first term in (39) yields the relative 
size of the effect of bending on the spring constant. A 
typical value of the relative change in spring constant 
due to bending was about 0.01 - 0.1% up to N = 200 

at sj\^r- = 0.05. At high N we would need a smaller 
deflection Z to reach the local maximum and minimum 
points in the curve in fig [21 however the relative change 
tends to be smaller as N is increased. From (39) and (35) 
we can calculate /3/(Cwq) in equation (35). This yields 



-4Nn 



y/N—n 



J__ 



(40) 



A plot of equation (40) is shown in fig [H 

What effect has a finite temperature on this result? 
Following Blom et alQ we have the grand canonical po- 
tential — E to t — (J-Ntot where the chemical potential 
fx w Ep at room temperature and 

/>oo 

Ntot = V / g(E-E n )f(E) dE, (41) 

POO 

Etot = V / g(E - E n )f(E)E dE, (42) 
z Je„ 
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FIG. 5: Example of spring constant as a function of different 
bending Z' = for four different temperatures using 

ekvation (43). Curve a: fc B T = 0.0001 eV, curve b: k B T = 
0.01 eV, curve c: feT = 0.025 eV (room temperature) and 
curve d: ksT = 0.05 eV. The plot is made using /i = 2.9 eV, 
Eq — 0.47 eV (i.e. do = 1 nm) and the upper energy limit of 
integration is taken to be E cut = 5.0 eV. Increasing E cut to 
7.0 yields no visible change in the curves. For E cut > n the 
contribution to the expression decreases rapidly due to the 
Fermi-Dirac function. 



where g(E) is the one dimensional density of states, f(E) 
is the Fermi-Dirac distribution function and E n is the 
energy of the state n that fits the cross section. Using 
F = — J§ = kZ and equation (32) we argue that the 
generalised expression for the spring constant k valid for 
any temperature should be given by 
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(44) 
(45) 



and the upper integration limit E cu t — > oo. Equation 
(43) reduces to (32) when T — > 0. A plot of equation 
(43) for different temperature is shown in fig [5] 

Due to temperature, more states becomes available for 
the electrons from E F up to about E F + k B T and de- 
creasing the number of available states between Ep and 
about Ep — ksT . This effect tend to make the transition 



at integer numbers in fig [4] (where a new energy level is 
added) less sharp, smoothing the spikes in £ = (3/(Cu)q). 
We also see in fig [5] that at room temperature the curve 
follows the zero temperature curve well, except around 
the minima where a new state is added due to bending. 
The room temperature curves here becomes smoothed. 
The condition that temperature are not important other 
than close to N being an integer must be that k B T <C Eq 
where Eq is the difference in energy between two energy 
levels. Using equation (26) we obtain this condition as 



knT knTrnd 2 , 



En 



2nh 2 



< 1. 



(46) 



To obtain agreement between the curve for finite tem- 
perature and the zero temperature curve around Z = 
we need equation (46) to be fullfilled. Room tempera- 
ture corresponds to = 5.3 % in fig [SJ This means 
that even for a diameter of 1 nm the system needs to 
be cooled-^ if one is to use the zero temperature result. 
However, interesting results can also be obtained at room 
temperature, as we see in the curves in fig[5j 

Experimentally, a way to detect these amplitude de- 
pendent resonance frequencies might be in situ TEM 
probing where one can see the wire while manipulat- 
ing it J^iH For the thermal vibration of a nanowire we 
have ifco^ 2 = kgT. For weak nanowires this ampli- 
tude is relatively large at room temperature and can be 
observed^ If one choose to drive the oscillation with 
an external electric field this effect must be taken into 
account. Small metallic nanosized cantilever has been 
manufactured^^— 



IV. CONCLUSIONS 

Using a one dimensional jellium model and standard 
beam theory we calculate the spring constant of a vi- 
brating nanowire cantilever. By using the asymptotic 
energy eigenvalues of the standing electron waves over 
the nanometer sized cross section area, the change in 
the grand canonical potential is calculated and hence the 
force and the spring constant. As the wire is bent more 
electron states fits in its cross section. This has an im- 
pact on the spring "constant" which oscillates slightly 
with the bending of the wire. In this way we obtain 
an amplitude dependent resonance frequency of the os- 
cillations that should be detectable. Because the weak 
non-linearity the resonance frequency is shifted from lvq 

to J loI + jf3A 2 . Using (40) we can replace j3 with £Cojq. 



We then obtain from this the relative frequency shift 



Aw _ 27 2 MA (A 



(47) 



where A is the amplitude of the oscillation and lj = 2irf. 
The data of some wires possible to use in an experiment 
is shown in table HI 
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TABLE I: Table of different L = 40 mn long gold nanowires (y = 0.44) that may be used to measure the predicted effects. The 
maximum temperature T of the wire and its diameter do are coupled by equation (46), if one wants to use the zero temperature 
result. However, as seen in fig [5] the modification of the curves due to temperature are rather small at room temperature and 
interesting measurements on the system can also be made at this higher temperature. £ varies periodically with increasing 
N — -g£ and is obtained from equation (40) for the zero temperature case. Small changes in do (yielding Eq) can result in 
large changes in £ if N is close to an integer. At about 25-30% of the distance between N being integers £ becomes zero as it 
change sign from positive to negative. The frequency shift is proportional to the square of the amplitude A. In the table we 
use A = 12 nm, that is 30% of the wires length L. For the thermal vibration of a nanowire we have ^koA 2 = ksT. Using this 
equation for the weakest nanowire in the table this amplitude becomes 13 nm at room temperature. 



d (nm) 


T (K) 


M(xl0 -22 kg) 


ko(mN/m) 


t 


fa (MHz) 


A/ (kHz) 


1.0 


55° 


1.93" 


0.0469 c 


-0.344" 


39. 2 e 


-0.992' 


1.5 


25" 


4.34' 


0.238 c 


-0.464" 


58.9 e 


-4.52 / 


2.0 


14° 


7.72 b 


0.750 c 


-0.395 d 


78.4 C 


-9.12 / 


2.5 


8.8° 


12. 1 6 


1.83 c 


-0.281"* 


97.8 C 


-12. & 


3.0 


6.2 a 


17.4 6 


3.80 c 


-0.107 d 


118 e 


-8.35 / 


3.5 


4.5 a 


23.6 6 


7.05 c 


-0.200 d 


138 e 


-24. 7 f 


4.0 


3.4 a 


30.1 b 


12. C 


-0.198" 


159 e 


-37.0 / 



"Prom equation (46) using fc J^ T =1% and do in the table. 

'The mass of the wire M = pLd 2 where p is the density. 
c Calculated using equation (34) and the data in the table. 
"Calculated using equation (40) and the data in the table. 



e From /o = using fco and M in the table. 

'From equation (47) using the data in the table. A = 12 nm. 
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